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A unified treatment of the cohesive and conducting properties of metalHc nanostructures in terms 
of the electronic scattering matrix is developed. A simple picture of metallic nanocohesion in which 
conductance channels act as delocalized chemical bonds is derived in the jellium approximation. 
Universal force oscillations of order ep /^f are predicted when a metallic quantum wire is stretched 
to the breaking point, which are synchronized with quantized jumps in the conductance. 

PACS numbers: 03.40.Dz, 62.20.Fe, 73.20.Dx, 73.40. Jn 

Cohesion in metals is due to the formation of bands, which arise from the overlap of atomic orbitals. In a metallic 
constriction with nanoscopic cross section, the transverse motion is quantized, leading to a finite number of subbands 
below the Fermi energy ep- A striking consequence of these discrete subbands is the phenomenon of conductance 
quantization The cohesion in a metallic nanoconstriction must also be provided by these discrete subbands, which 
may be thought of as chemical bonds which are delocalized over the cross section. In this Letter, we confirm this 
intuitive picture of metallic nanocohesion using a simple jellium model. Universal force oscillations of order ep/^F 
are predicted in metallic nanostructures exhibiting conductance quantization, where Ai? is the Fermi wavelength. Our 
results are in quantitative agreement with the recent pioneering experiment of Rubio, Agrai't, and Vieira who 
measured simultaneously the force and conductance during the formation and rupture of an atomic-scale Au contact. 
Similar experimental results have been obtained independently by Stalder and Diirig |^ . 

Quantum-size effects on the mechanical properties of metallic systems have previously been observed in ultrasmall 
metal clusters Q, which exhibit enhanced stability for certain magic numbers of atoms. These magic numbers have 
been rather well explained in terms of a shell model based on the jellium approximation Q. The success of the 
jellium approximation in these closed nanoscopic systems motivates its application to open (infinite) systems, which 
are the subject of interest here. We investigate the conducting and mechanical properties of a nanoscopic constriction 
connecting two macroscopic metallic reservoirs. The natural framework in which to investigate such an open system 
is the scattering approach developed by Landauer j|] and Biittiker |6|. Here, we extend the formalism of Ref. 
which describes electrical conduction, to describe the mechanical properties of a confined electron gas as well. 

For definiteness, we consider a constriction of length L in an infinitely long cylindrical wire of radius R, as shown 
in Fig. |l|. We neglect electron-electron interactions, and assume the electrons to be confined along the z axis by a 
hard- wall potential at r = r(z). This model is considerably simpler than a self-consistent jellium calculation but 
should suffice to capture the essential physics of the problem. Outside the constriction, the Schrodinger equation is 
separable, and the scattering states can be written as 

^LJ'/'- ^) = e±^'=-+™* J,„(7™„r/i?), (1) 

where the quantum numbers 7mn are the roots of the Bessel functions Jmilmn) = 0. These scattering states may be 
grouped into subbands characterized by the quantum numbers m and n, and we shall use the notation v = (m, n). 
The energy of an electron in subband v is £{k) = 6,^ + fi^k'^ /2m, where 



(2) 



The fundamental theoretical quantity is the scattering matrix of the constriction S{E), which connects the incoming 
and outgoing scattering states. For a two-terminal device, such as that shown in Fig. |l|, S{E) can be decomposed into 
four submatrices Sap{E), a, P = 1,2, where 1 (2) indicates scattering states to the left (right) of the constriction. 
Each submatrix SapiE) is a matrix in the scattering channels vv' . 

In terms of the scattering matrix, the electrical conductance is given by 

G=^J dE^Tr{sUE)S,2iE)}, (3) 

where f{E) = {exp[/3(i? — fi)] + is the Fermi distribution function and a factor of 2 has been included to account 
for spin degeneracy. The grand canonical potential of the system is 
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n = -keT J dE D{E) In [l + e-'^^^-i'^'j , (4) 
where the density of states in the constriction may be expressed in terms of the scattering matrix as 0^ 

Eqs. (|^) to allow one to treat the conducting and mechanical properties of a confined electron gas on an equal 
footing, and provide the starting point for our calculation. 

We are interested in the mechanical properties of a metallic nanoconstriction in the regime of conductance quan- 
tization. The necessary condition to have well-defined conductance plateaus in a three-dimensional constriction was 
shown by Torres, Pascual, and Saenz to be {dr/dz)"^ <C 1. In this limit, Eqs. (^) to (||) simplify considerably because 
one may employ the adiabatic approximation u^. In the adiabatic limit, the transverse motion is separable from the 
motion parallel to the z axis, so Eqs. (0) and (g) remain valid in the region of the constriction, with R replaced by 
r(z). The channel energies thus become functions of z, £,j{z) = /i^7^/2mr(z)^. In this limit, the scattering matrices 
SapiE)^ a, /3 = 1, 2 are diagonal in the channel indices, leading to an effective one-dimcnsional scattering problem. 
The condition (dr jdzY 1 and the requirement that the radius of the wire outside the constriction not be smaller 
than an atomic radius (i.e., kpTl > 1) automatically imply the validity of the WKB approximation. 

Since the energy differences between the transverse channels in an atomic-scale constriction are large compared 
to fcsT at ambient temperature, we restrict consideration in the following to the case T = 0. In the adiabatic 
approximation, the conductance becomes 

where the transmission probability for channel v may be calculated using a variant of the WKB approximation ||ll|,|l2| , 
which correctly describes the rounding of the conductance steps at threshold. The density of states in the constriction 
in the adiabatic approximation is 

V 

where the total phase shift is given in the WKB approximation by 

e,{E) = {2m/n''y^ / dz[E- e,{z)]^/^ (8) 

"'0 

the integral being restricted to the region where e^{z) < E. The grand canonical potential of the system is thus 

the sum being over channels with e,y{z) < ep- Under elongation, the tensile force is given by F = —dVt/dL. It is easy 
to show that F is invariant under a stretching of the geometry r[z) ~f r(Az), i.e., 

F^^f{AL/Lo,kFR), (10) 

where f{x,y) is a dimensionless function. Nonuniversal corrections to F occur in very short constrictions, for 
which the adiabatic approximation breaks down. The leading order correction to the integrand in Eq. (^) is 
— (37r/64)fci?r(z)((ir/(iz)^, leading to a relative error in F of ~ 2sin'^6'/4, where 9 is the opening angle of the con- 
striction. Using a modified Sharvin equation |^ to estimate the diameter of the contact versus elongation for the 
experiment of Ref. ||^ indicates an opening angle 9 < 45°, for which the nonuniversal corrections are ^ 8%, justifying 
the above approach. 

Fig. H shows the conductance and force of a metallic nanoconstriction as a function of the elongation, calculated 
from Eqs. (^) and (|^). Here an ideal plastic deformation was assumed, i.e., the volume of the constriction was 
held constant [|l3|. The correlations between the force and the conductance are striking: \F\ increases along the 



2 



conductance plateaus, and decreases sharply when the conductance drops. The constriction becomes unstable when 
the last conductance channel is cut off. Some transverse channels are quite closely spaced, and in these cases, the 
individual conductance plateaus [e.g., G/(2e^//i) = 14, 15, 19, 21] and force oscillations are difficult to resolve. Fig. || 
is remarkably similar to the experimental results of Refs. |^ and both qualitatively and quantitatively. Inserting 
the value Ef/^f — 1.7nN for Au, we see that both the overall scale of the force for a given value of the conductance 
and the heights of the last two force oscillations are in quantitative agreement with the data shown in Fig. 1 of Rcf. 

We wish to emphasize that the calculation of F presented in our Fig. || contains no adjustable parameters [ p4[ . 
The increase of |F| along the conductance plateaus and the rapid decrease at the conductance steps were described 
in Ref. Q as "elastic" and "yielding" stages, respectively. With our intuitive picture of a conductance channel as a 
delocalized metallic bond, it is natural to interpret these elastic and yielding stages as the stretching and breaking of 
these bonds. 

The fluctuations in F due to the discrete transverse channels may be thought of as arising from finite-size corrections 



to the surface tension a. However, as in the case of universal conductance fluctuations |15|, it is more instructive to 
consider the extensive quantity F itself, rather than the intensive quantity cr. Approximating the sum in Eq. (^) by 
an integral, and keeping the leading order corrections, one obtains 

n = ujV + c7S-—^L + sn, (11) 

3A F 

where V is the volume of the system, S is the surface area, lu = -~2eFkp/15Tr'^ is the macroscopic free energy density, 
and a = ei?fc|,/167r is the macroscopic surface energy. The remaining term SQ is a quantum correction due to the 
discrete transverse channels, and may be either positive or negative. Under an ideal plastic deformation, the volume 
of the system is unchanged, and the tensile force is 

where SF — —d{5^)/dL. The ffi'st term in Eq. dl3) is the contribution to the force due to the macroscopic surface 
tension. This is plotted as a dashed line in Fig. l2l for comparison. The macroscopic surface tension determines the 
overall slope of F . The quantum corrections to F due to the discrete transverse channels consist of a constant term 
plus the fluctuating term 5F. Fig. |^ shows 5F for three different geometries and for values of RfR from 6 to 1200, 
plotted versus the corrected Sharvin conductance ||^ 

^ — , (13j 

where Amin and Cmin are the area and circumference of the constriction at its narrowest point. Gs gives a smooth 
approximation to G. As shown in Fig. ^(a), the force oscillations obey the approximate scaling relation 

6F{AL/Lo, kFR) ^ ^Y(G.), (14) 

Af 

where y is a dimensionless scaling function which is independent of the precise geometry r{z). Eq. ( [l^ ) indicates 
that the force fluctuations, like the conductance, are dominated by the contribution from the narrowest part of the 



constriction, of radius i?min. The scaling relation (14) breaks down when Rmin/R ^ 0.8. 

Fig. H shows that the amplitude of the force fluctuations persists essentially unchanged to very large values of Gs 
It was found to be 

— 2\ 1/2 



[Y^-Yj -0.3 (15) 



for < Gg < 10^. The detailed functional form of Y{Gs), like the distribution of widths of the conductance plateaus, 
depends on the sequence of quantum numbers 7^ , which is determined by the shape of the cross section. However, the 
amplitude of the force fluctuations AY was found to be the same for both circular and square cross sections. Both 
these geometries are integrable, and hence have Poissonian distributions of transverse modes. It is clearly of interest 
to investigate the force fluctuations for nonintegrable cross sections, with non-Poissonian level statistics. 

The experiments of Refs. [|[|| observed well-defined conductance steps, but found no clear evidence of conductance 
quantization for G/(2e^//i) > 4 ||l^. Deviations of the conductance plateaus from integer values in metallic point 
contacts are likely to be due to backscattering from imperfections in the lattice or irregularities in the shape of the 
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constriction . We find that such disorder- induced coherent backscattering leads to noise-hke fine structure in 
the conductance steps and force oscillations, with a reduction of the conductance on the plateaus, but no shift of the 
overall force oscillations Our prediction of universal force oscillations is consistent with the experiments of Refs. 
1^ and H], which found force oscillations with an amplitude comparable to our theoretical prediction for G/{2e^ /h) 
up to 60. 

Molecular dynamics simulations by Landman et al. Todorov and Sutton |20|, and Brandbyge et al. have 
suggested that the conductance steps and force oscillations observed in Refs. and |^ may be due to a sequence 
of abrupt atomic rearrangements. While the discreteness of the ionic background is not included in the jellium 
model, our results nevertheless suggest that such atomic rearrangements may be caused by the breaking of the 
extended metallic bonds formed by each conductance channel. However, it should be emphasized that our prediction 
of universal force fluctuations of order ef/^f is not consistent with the simulations of Refs. and [^, which predict 
force fluctuations which increase with increasing contact area. This discrepancy may arise because we consider the 
equilibrium deformation of a system with extended electronic wavefunctions, while Refs. |12 19|j20[] use a purely local 
interatomic potential and a fast, nonequilibrium deformation 

In conclusion, we have presented a simple jellium model of metallic nanocohesion in which conductance channels act 
as delocalized metallic bonds. This model predicts universal force oscillations of order ef/^f in metallic nanostructurcs 
in the regime of conductance quantization, and is able to explain quantitatively recent experiments on the mechanical 
properties of nanoscopic metallic contacts H,^. The formalism developed here based on the electronic scattering 
matrix should be applicable to a wide variety of problems in the rapidly evolving fleld of nanomechanics. 
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FIG. 1. Schematic diagram of a constriction in a cylindrical quantum wire. Electrons are confined along the z axis by a 

hard-wall potential at r = "riz)- Two different geometries are considered: r(z) = {R + _Rmin)/2 + {R — Rmin) cos(27rz/L)/2 
(cosine constriction) and r{z) = Rmin + (R — Rmin){2z / L — Vf (parabolic constriction), with r{z) = R ior z < Q and z > L. 
The minimum radius of the neck i?min as a function of the elongation AL/Lo is determined by a constant volume constraint 
J^\{zrdz = R'Lo. 
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FIG. 2. Electrical conductance G and tensile force F of a cosine constriction in a cylindrical quantum wire of radius 
kpR = 11 versus the elongation AL/Lo- For the calculation of G, an initial length kpLo = 50 was assumed. The dashed line 

indicates the contribution to the force due to the macroscopic surface tension Fs ~ —adS/dL, whore S is the surface area of 
the system and a = epk^/Wn. Fs determines the overall slope of F, on which are superimposed the quantum oscillations due 
to the discrete transverse channels. 
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FIG. 3. Force fluctuations versus the corrected Sharvin conductance, (a) Results for cosine and parabolic constrictions in 
cylindrical quantum wires with radii kpR = 6, 8, 10, 11, and 12 (left to right). The results for the cosine and parabolic 

constrictions arc almost indistinguishable when plotted as a function of Gs , despite the fact that the total elongation and total 
work done differ by roughly a factor of 2 in the two cases. The dependence on R is also very weak, except for J?min « R (the 
rightmost portion of each curve), (b) Force fluctuations for a cosine constriction in a square quantum wire with 120 conductance 
channels, and (c) for a wire with 1200 conductance channels. 
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